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ABSTRACT 

Strange quark stars (SSs) may originate from accreting neutron stars (NSs) in low- 
mass X-ray binaries (LMXBs). Assuming that conversion of NS matter to SSs occurs 
when the core density of accreting NS reaches to the density of quark deconfinement, 
~ 5p0) where po ~ 2.7 x 10 14 g cm~ 3 is nuclear saturation density, we investigate LMXBs 
with SSs (qLMXBs). In our simulations, about l%o — 10% of LMXBs can produce SSs, 
which greatly depends on the masses of nascent NSs and the fraction of transferred 
matter accreted by the NSs. If the conversion does not affect binaries systems, LMXBs 
evolve into qLMXBs. We find that some observational properties (spin periods, X-ray 
luminosities and orbital periods) of qLMXBs are similar with those of LMXBs, and it 
is difficult to differ them. If the conversion disturbs the binaries systems, LMXBs can 
produce isolated SSs. These isolated SSs could be submillisecond pulsars, and their 
birthrate in the Galaxy is ~5-70 per Myr. 

Subject headings: b 
inaries: close — stars: neutron — dense matter 



1. Introduction 



There are at least three different ki nds of c o mpac t stars in the Universe: white dwarfs (WDs), 
neutron stars (NSs), and black holes. IWittenl ()1984l ) suggested a possible existence of compact 
objects consis t ing o f str ange quark matter. Due to strange quark matter being absolutely stable, 



Haensel et al.l (|1986l ) and lAlcock et a 



matter and not neutrons. However, 



( 1986 ) poin ted out that NSs almost would be made of strange 



Alpar (1987) considered that glitching radio-pulsars are NSs 



and not strange quark stars (SSs). iMadsenl (jl988l ) suggested that SSs can not b e formed directly 
in supernovaqj, or less they would eventually contaminate the entire Galaxy. iKluzniakl ( 1994 ) 
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Dai et al 



( 1995T ) and IXu et al.l l|200lf ) suggested that SSs can be formed directly during or shortly after some 



supernovae explosion when the central density of the proto-NSs is high enough. 



- 2 - 



suggest ed that SSs cou ld exist as millisecond pulsars. Due to the fast rotation and thermonuclear 
bursts, iLi et al.l (jl999h suggested that the SAX J1808.4-3658 is a good SS candidate. These SSs 
can be formed in low-ma ss X-ray binaries (L MXBs) via an accretion-triggered phase transition of 
NS matter to SS matter (jcheng & Dailll996l ). 



The phase transition requires the formation of a strange matte r seed in the NS. Th e strange 
matter is produced through the neutron matter at a critical density. ISerot &: Uechil (|1987l ) pointed 
out that the central density of an 1.4 M NS with a rather stiff equation of state is sufficientl y 
lower than the criti cal density. Based on the modern equations of state in IWiringa et al.l (|1988), 
Cheng &: Pail (| 19961 ) estimated that the NSs with 1.4 M Q must accrete matter of ~ 0.5M Q in order 
that their central densities reach the deconfinement density. Once the above condition is satisfied, 
the phase transition occurs. 



Olintol (|1987l ) prop osed that the process of the st range matter swallowing the neutron matter is 
a slow mode. Howev er, iHorvath fe; Benvenutd (|1988l ) showed that it is hydrodynamically unstable. 
Cheng Pail (| 19961 ) proposed that the conversion of neut ron matter sh o uld p roceed in a detonation 
mode and could be accompanied by a gamma-ray burst. lOuved et al.l (|2002l ) suggested that there 
is a quark-nova when the core of a NS (having experienced a transition to an up and down quark 
phase) shrinks into the equilibrated quark object after reaching strange quark matter saturation 
density (where a composition of up, down and strange quarks is the favored state of matt e r). In 



their model, the energy released as radiation in a quark-nova is up to 10 53 ergs. lOuved et al.1 (120111 



proposed that the quark novae in LMXBs may be the engines of short gamma-ray bursts. 

Based on the above descriptions, it is possible that SSs originate from the hydrodynami- 
cally unstable conversio n or the sl o w con version in LMXBs. Using standard equation of states 
of neutron-rich matter, IStaff et al.l (|2006l ) considered that the density of quark deconfinement is 
~ 5po, w here pn ~ 2.7 x 10 1 4 g cm -3 is nuclear saturation density. According to the equation of 
states in 



Akmal et al 



[4998), the gravitational mass of a NS is ~ 1.8M@ in order to reach 5po- 
Therefore, it is very important for our understanding of SSs' formation to study the mass evolution 
of NSs in LMXBs. 

In this work, by simulating the interaction of a magnetized NS with its environment and 
utilizing a population synthesis code, we focus on the mass change of NSs in LMXBs and the 
possibility from NSs converting SSs in LMXBs, and investigate the properties of LMXBs with SSs 
(qLMXBs). In Section 2, we present our assumptions and describe some details of the modelling 
algorithm. In Section 3, we discuss the main results and the effects of different parameters. In 
Section 4, the main conclusions are given. 



2. Model 



For the simulation of bin ary evolution, we use r apid binary star evolution code BSE (jHurlev et al 
2002) which was updated bv lKiel fc Hurlevl (|2006l ). In interacting binaries, NSs can be formed via 
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three channelsfe.g.. Ilvanova et alj|200a : Kiel et al.ll2008l ): (i) Core-collapse supernovae (CCSN) for 
a star; (ii)Evolution induced collapse (EIC) of a helium star with a mass between 1.4 and 2 .5M^ 
in which the collapse is triggered by electron capture on 20 Ne and 24 Mg ( Mivaii et al.lll980l ): (hi) 
Accretion-induced collapses (AIC) for an accreting ONeMg WD whose mass reaches the Chan- 
drasekhar limit. Respo nse of accreting ONeMg WD is treated in the same way as the evolution of 
CO WD (see details in iLu et alJlionflf ). 



2.1. Mass of Nascent NS 



Possibly the mass is one of the most important properties of NSs. However, the mass distribu- 
tion of nascent NSs is not yet well known. In BSE code, the gravitati onal mass of a nasc ent NS via 
CCSN depends on the mass of the CO-core at the time of supernova (jHurlev et al.ll2000l ). Figure Q] 
shows the masses of nascent NSs forming from different ini tial masses. Some authors assumed that 
the initial masses of NSs (Mj^ Tg 



2002 



are 1.4 Mq in their works (Ergm a et al 



Nelson Rappaportll200a . e.g.,). Lattimer Sz Prak ash (2007) showed that the masses of some 



1998 



Podsiadlowski et al 



NSs are lower than 1.4 Mq. Recently, Ivan der Meer et al.l (J2007J) found that the masses of NSs in 
SMC X-l and Cen X-3 are 1.061^^ M © and L34 -ai4 M ©' respectively. However, It is well known 
that most of the accurately measured masses of NSs are near 1.4 Mq. 



In our work, we use the initial masses of NSs via CCSN in iHurlev et al.l (120000 and M l^ = 
1.4Mq in different simulations, respectively. For NSs via AIC, following IHurlev et al.l ((2000), we 



take M^ s 



1.3M . Similarly, for NSs via EIC, we also take ML = 1.3M, 



i . 



In addition, nascent NS receives additional velocity ("kick") due to some still unclear process 
that disrupts spherical sym metry during the collapse or later Dichotomous nature of kicks which 
was suggested quite early by iKata ()1975l ) . Observationally, the kick is not well constrained due to 
numerous selection effects. Currently, high kicks (~ 100 km s _1 ) are associa ted with NS origin ating 
from CCSN, while low kicks (~ 10km s" 1 ) with NS born in EIC and AIC (jpfahl et alJbood ). We 
apply to core-collapse NS Maxwellian distribution of kick velocity 



2vl 
V 



(1) 



Variation of a k between 50 and 200 k m s -1 , introduces a n uncertain t y < 3 in the birthrate of low- 
and intermediate- mass X-ray binaries (IPfahl et al.ll2003l ). IZhu et al.1 (|2012l ) discussed the effects of 
parameter o"k on LMXBs' populations. Since in this paper we focus on the physical parameters that 
mostly a f fect th e masses of NSs, we do not discuss the effects of on SSs' population. Following 
Lii et al.l (j2012h . we take a k = 190 km s" 1 in CCSN, and cr k = 20 km s' 1 in EIC and AIC. 
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2.2. Mass of Accreting NS 



In LMXBs, NSs accrete the matter from their companions via Roche lobe flows or stellar winds. 
The interaction of a rotating mag netized NS (single or in a binary system) with surrounding matter 
has been studies by many author s ( Pringle Sz Reesfl972 ; Illarionov Sunvaevil975 ; Ghosh &: Lamb 
19781 : lLovelace et al.lllQQfj . Il99fll . e. g.,). 

Using a convenient way of describing NS evolution elaborated by Lipunov et al. J 19921) and 



a. 



recent model for quasi- spherical accretion including subsonic settling proposed by iShakura et al 
(J2012J), iLii et al.l (|2012l ) gave detailed simulations for sp in period evolut ion and matter accretion 
of NSs in binaries. In this work, we adopt their model. ILii et al.l (|2012l ) assumed that all matter 
transferred is accreted by the NS in an accretor stage. We introduce a parameter f3 which is the 
fraction of transferred matter accreted by the NS, and the rest of the transferred matter is lost from 
binary system. The lost matter takes away t he specific angular momentum of t he prospective donor. 



1992 



The value of /3 has been usually set to 0.5 (jPodsiadlowski et al. 
2003). In our work, we set (3 = 0.1, 0.5 and 1.0 in different simulations 



20021 ; iNelson fc Rappaport 



2.3. LMXBs with Nascent SSs 



As noted in Jfl when the gravitational mass of a NS reaches to ~ 1.8Mq, the NS can turn into 
a SS via the hydrodynamically un stable conversion or t he slow conversion. In the former mode, the 
binary system may be disrupted ()Ouved fc Stafljboill ). and it becomes two isolated stars. One of 



them is an isolated SS. However, it is difficult to know the effects of the hydrodynamically unstable 
conversion on binary systems. Therefore, in our work, we consider two extreme cases: (i) In order 
to simulate all potential properties of qLMXBs, we assumed all LMXBs are not affected and survive 
after the conversion, that is, LMXBs become qLMXBs when the masses of accreting NSs are larger 
than 1.8M ; (ii) We assumed that all LMXBs are disrupted after the conversion, that is, there is 
no qLMXBs. However, we can discuss the origin of isolated submillisecond pulsars. 

In our simulations, if a nascent NS has larger mass than 1.8 Mq it is a SS. Therefore, in the 
paper, qLMXBs include some LMXBs in which the nascent NSs have larger masses tha n 1.8 M^. 



Lai Xul (120091 ) suggested that SSs could have high maximum masses (See Figure 2 in lLai Xul 
(|2009l )) if they are composed of the Lennard-Jones matter. In our work, we assume that the 
maximum mass of SS is 3.0 Mq. 



3. Results 



We use Monte Carlo method to simulate the initial binaries. For initial mass function, mass- 
ratios, and separation s of components in binary systems, we adopt the distributions used by us in 
Lii et aJj (|2006l . l2008l ). We assume that all binaries have initially circular orbits. After a supernova, 
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new parameters of the orbit are derived using standard formulae(e. g., iHurlev et al.ll2002l ). It is 



well known that theoretical models of the population of LMXBs depen d on badly known input pa- 



rame ters, such as kick velocity and common envelope treatment (e. g., iPfahl et al.ll2003l ; IZhu et al 



2012 ) . However, in this pioneering study of qLMXBs we focus on the effects which are important 
for the masses of NSs: the efficiency accreted of transferred matter, /? ( f3 = 0.1, 0.5 and 1.0 ), 
and the mass of nascent NS via CCSN (See Figure [T] ). We use 1 x 10 8 binary systems in our 
Monte-Carlo simulations. 



3.1. Mass Increases of Accreting NSs 

Figure [2] gives the mass increases of accreting NSs in LMXBs. According the assumptions of 
SSs formatting from NSs with masses higher than 1.8 M , l%o {fi = 0.1) — 10% (f3 = 1.0) of 
LMXBs are qLMXBs in our simulations. This proportion greatly depends on input parameters /3 
and the initial masses of nascent NSs. In the cases of j3 = 1.0 and /3 = 0.5, most of SSs in qLMXBs 
come from low-mass NSs with 1.4 M , that is, they have accreted ~ O.4M matter. In the case of 
f3 = 0.1, most of SSs in qLMXBs originate from CCSN. 

Mass increases of accreting NSs in LMXB s depend on not only input parameter /3 in this work, 



but also orbital periods and NSs' companions. IZhu et al.l (|2012l ) showed that most of NSs in LMXBs 



with WD donors have low mass-accretion rates (~ 10 12 M yr ) and most of LMXBs with WD 
donors are transient. Therefore, the masses of accreting NSs in LMXBs with WD donors hardly 
reach to 1.8 M Q . Less than 10% (in case of /3 = 0.1)— l%c (in case of /3 = 0.5, M^ s = 1.4M ) of 
qLMXBs have WD donors in our simulations. Most of qLMXBs have main sequence donors. From 
now on, we just discuss LMXBs and qLMXBs with MS donors. 



3.2. Properties of qLMXBs 

According to our assumption that the conversion of NS matter to SS matter does not affect 
binary systems, we can simulate some observational properties of LMXBs and qLMXBs, and wish 
to find some differences between them. Then, taking the case of /3 = 0.5 as an example, we discuss 
some properties of qLMXBs. 

X-ray luminosities (mass-accretion rates), spin periods and orbital periods are important pa- 
rameters of LMXBs. Figure [3] shows the accretion rates by NSs or SSs in LMXBs and qLMXBs, 
and the X-ray luminosities which are approximated as 

L x = t/Mnsc 2 ~ 5.7 x 10 35 erg s" 1 ^) x ( 10 -ife?| yr -i ), (2) 

where rj ~ 0.1 is the efficiency of converting accreted mass into X-ray photons. We can find that 
there is not significant difference between accretion rates by NSs or SSs in LMXBs and qLMXBs. 
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Fig. 1. — Masses o f nascent NSs via CCSN vs. stellar initial masses. The solid line comes from 
Hurley et al. (j200ol ). and the dashed line means that masses of nascent NSs are equal to IAMq. 
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Fig. 2. — Distributions of the final masses of compact stars (NSs or SSs) vs. their initial masses in 
LMXBs for different input parameters. Gradations of gray-scale correspond to the number density 
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The spin p eriods of the rotating magnetized NSs mainly depend on their mass-accretion rates 
(|Lii et al.ll2012l ). If accreting SSs are similar to NSs, we can simulate the spin periods of accreting 
SSs. Figure S] gives the distribution of the spin periods of NSs in LMXBs or SSs in qLMXBs. In 
general, spin periods of SSs in qLMXBs are shorter than those of NSs in LMXBs. 

Figure [5] shows the the distributions of the orbital periods P Q rb of qLMXBs and LMXBs vs. 
masses of their secondary stars. Similarly, there is not significant difference between qLMXBs and 
LMXBs. 

As the above descriptions show, it is very difficult in our model to differ from LMXBs and 
qLMXBs. The most effective way is to measure the masses of compact objects in LMXBs if the 
assumption that NSs with masses larger than 1.8 Mq are SSs is right. 



Jonker et al.l (|2005l ) suggested that compact object in 2S 0921-630 (It is a LMXB) has a 
mass between ~1.9 — 2.9 Mq. Accordin g to our assumption , the compact object is a SS. However, 



Steeghs &: Jonkerl ([20071 ) considered that iJonker et al.l ([2009 ) overestimated the rotational broaden- 
ing and the mass of compact object in 2S 0921-630 is ~1.44 M & . The orbital period of 2S 0921-630 



is 9.006±0.007 days, and its X-ray luminosity is ~ 10 3b erg s _1 ([Kallman et al J 1200.1 ). Results of 



simulating LMXBs and qLMXBs in our work both cover these observations. Therefore, we can not 
conclude whether compact object in 2S 0921-630 is NS or SS. 



Pernor est et al.l (|2010l ) gave that PSR 1614-2230 has a mass of 1.97 Mq. And, it is a millisecond 



radio pulsar (Pulsar spin period is 3.15 ms) in an 8.7 day orbit, and its companion has a mass of 



0.5 M q. Although PSR 1614-2230 is not X-ray binary, it may come from X-ray binary. iLin et al 



([201 lh suggested that PSR 1614-2230 desce nded from a LMXB v ery much like Cyg X-2 (P orb =9.8 



days, Mns = 1.7M© and M2 = O.6M , see ICasares et al.l (I2010l )). AS Figure [5] shows, our results 



cover the positions of orbital periods and companion masses of Cyg X-2 and radio millisecond 
pulsar binary PSR 1614-2230. Our work support that PSR 1614-2230 originate from a LMXB. 
PSR 1614-2230 may be a SS. 



3.3. Submillisecond Pulsars 



Weberl (|2005l ) suggested that an isolated submillisecond pulsar spinning at ~ 0.5 ms could 
strongly hint the existence of SSs. If the conversion of NS matter to SS matter is hydrodynamically 
unstable and LMXBs are disrupted, the nascent S Ss are isolated. Assuming; that one binary with 



Mi > 0.8M Q is formed per year in the Galaxy (jYungelson et al 



1993 



Han 



1998 



Hurley et al 



20021 ). we can estimate that the occurrence rate of hydrodynamically unstable conversion is about 
5—70 per Myr. 

As Figure [6] shows, the majority of NSs at the beginning of the conversion have spin periods 
longer than 1 ms. If the angular momentum is conserved and no unknown physical mechanic 
spins up NSs and nascent SSs during the conversion, the spin periods of nascent SSs depend on 
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Fig. 3. — Distributions of accretion rates (X-ray luminosities) by NSs or SSs in LMXBs and 
qLMXBs, respectively. 
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Fig. 4. — Number distributions of the spin periods of NSs or SSs in LMXBs and qLMXBs, respec- 
tively. The numbers are normalized to 1. 
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the change of moment of inertia. lOuved &: Staff! (|201ll ) considered that a typical ejected mass of 
hydrodynamically unstable conversion is ~ 1O~ 3 M0. Therefo re, the chang e of m oment of inertia 
is determined by the difference of radius between NS and SS. lOuved et al.l (j2002l ) estimated that 
NS could shrink by as much as 30%. Then, many nascent SSs have spin periods of ~ 0.5 ms in our 
simulations. However, this result greatly depends on the equation of state of NSs and SSs which 
are poorly known. If NS only shrinks by as much as 10% (Private discussion with Xu), it is difficult 
for the nascent SSs to spin up to 0.5 ms. 



4. Conclusions 

Employing the population synthesis approach to the evolution of binaries and using the inter- 
acting model of a rotating magnetized NS with surrounding matter, we investigate the mass change 
of NSs in LMXBs and the possibility from NSs converting SSs in LMXBs. Our results show that 
about l%o — 10% of LMXBs can produce SSs. These SSs may exist in qLMXBs or be isolated, 
which depends on physical model of the conversion of NS matter to SS mater. 

Our toy model can not conclude whether there are SSs in the Galaxy and can not give what 
properties qLMXBs have. In further work, we need detailed physical model (equation of state 
about NS and SS and the conversion of NS matter to SS mater) to improve our work. 
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Low-mass X-ray Binaries with Strange Quark Stars 



Chunhua Zhu 1 ' 2,3 , Guoliang Lii 1,2 , Zhaojun Wang 1 ' 2 , Jinzhong Liu 3 
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ABSTRACT 

Strange quark stars (SSs) may originate from accreting neutron stars (NSs) in low- 
mass X-ray binaries (LMXBs). Assuming that conversion of NS matter to SSs occurs 
when the core density of accreting NS reaches to the density of quark deconfinement, 
~ 5po> where po ~ 2.7 x 10 14 g cm -3 is nuclear saturation density, we investigate LMXBs 
with SSs (qLMXBs). In our simulations, about l%o — 10% of LMXBs can produce SSs, 
which greatly depends on the masses of nascent NSs and the fraction of transferred 
matter accreted by the NSs. If the conversion does not affect binaries systems, LMXBs 
evolve into qLMXBs. We find that some observational properties (spin periods, X-ray 
luminosities and orbital periods) of qLMXBs are similar with those of LMXBs, and it 
is difficult to differ them. If the conversion disturbs the binaries systems, LMXBs can 
produce isolated SSs. These isolated SSs could be submillisecond pulsars, and their 
birthrate in the Galaxy is ~5-70 per Myr. 

Subject headings: b 
inaries: close — stars: neutron — dense matter 

1. Introduction 

There are at least three different kinds of compact stars in the Universe: white dwarfs (WDs), 
neutron stars (NSs), and black holes. ? suggested a possible existence of compact objects consisting 
of strange quark matter. Due to strange quark matter being absolutely stable, ? and ? pointed out 
that NSs almost would be made of strange matter and not neutrons. However, ? considered that 
glitching radio-pulsars are NSs and not strange quark stars (SSs). ? suggested that SSs can not 
be formed directly in supernova^], or less they would eventually contaminate the entire Galaxy. ? 
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? and ? suggested that SSs can be formed directly during or shortly after some supernovae explosion when the 
central density of the proto-NSs is high enough. 
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suggested that SSs could exist as millisecond pulsars. Due to the fast rotation and thermonuclear 
bursts, ? suggested that the SAX J1808. 4-3658 is a good SS candidate. These SSs can be formed 
in low-mass X-ray binaries (LMXBs) via an accretion-triggered phase transition of NS matter to 
SS matter (?). 

The phase transition requires the formation of a strange matter seed in the NS. The strange 
matter is produced through the neutron matter at a critical density. ? pointed out that the central 
density of a 1.4 M NS with a rather stiff equation of state is sufficiently lower than the critical 
density. Based on the modern equations of state in ?, ? estimated that the NSs with a 1.4 M Q must 
accrete matter of ~ 0.5M Q in order that their central densities reach the deconfinement density. 
Once the above condition is satisfied, the phase transition occurs. 

? proposed that the process of the strange matter swallowing the neutron matter is a slow 
mode. However, ? showed that it is hydro dynamically unstable. ? proposed that the conversion of 
neutron matter should proceed in a detonation mode and could be accompanied by a gamma-ray 
burst. ? suggested that there is a quark-nova when the core of a neutron star (having experienced a 
transition to an up and down quark phase) shrinks into the equilibrated quark object after reaching 
strange quark matter saturation density (where a composition of up, down and strange quarks is 
the favored state of matter) . In their model, the energy released as radiation in a quark- nova is up 
to 10 53 ergs. ? proposed that the quark novae in LMXBs may be the engines of short gamma-ray 
bursts. 

Based on the above descriptions, it is possible that SSs originate from the hydrodynamically 
unstable conversion or the slow conversion in LMXBs. Using standard equation of states of neutron- 
rich matter, ? considered that the density of quark deconfinement is ~ 5po, where po ~ 2.7 x 10 14 g 
cm~ 3 is nuclear saturation density. According to the equation of states in ?, the gravitational mass 
of a NS is ~ 1.8M Q in order to reach 5po- Therefore, it is very important for our understanding of 
SSs' formation to study the mass evolution of NSs in LMXBs. 

In this work, by simulating the interaction of a magnetized NS with its environment and 
utilizing a population synthesis code, we focus on the mass change of NSs in LMXBs and the 
possibility from NSs converting SSs in LMXBs, and investigate the properties of LMXBs with SSs 
(qLMXBs). In Section 2, we present our assumptions and describe some details of the modelling 
algorithm. In Section 3, we discuss the main results and the effects of different parameters. In 
Section 4, the main conclusions are given. 

2. Model 

For the simulation of binary evolution, we use rapid binary star evolution code BSE (?) with 
updates by ?. In interacting binaries, NSs can be formed via three channels(e.g., ??): (i) Core- 
collapse supernovae (CCSN) for a star; (ii)Evolution induced collapse (EIC) of a helium star with 
a mass between 1.4 and 2.5M© in which the collapse is triggered by electron capture on 20 Ne and 
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24 Mg (?); (iii) Accretion-induced collapses (AIC) for an accreting ONeMg WD whose mass reaches 
the Chandrasekhar limit. Response of accreting ONeMg WD is treated in the same way as the 
evolution of CO WD (see details in ?). 

2.1. Mass of Nascent Neutron Star 

Possibly the masses are one of the most important properties of NSs. However, the mass 
distribution of nascent NSs are not yet well known. In BSE code, the gravitational mass of a 
nascent NS via CCSN depends on the mass of the CO-core at the time of supernova (?). Figure [Q 
shows the masses of nascent NSs forming from different initial masses. Some authors assumed that 
the initial masses of NSs (M^ s ) are 1.4 Mq in their works (???, e.g.,). ? showed that the masses 
of some NSs are lower than 1.4 Mq. Recently, ? found that the masses of NSs in SMC X-l and 
Cen X-3 are 1.06to'iJ Mq and 1.34^ 014 Mq, respectively. However, It is well known that most of 
the accurately measured masses of NSs are near 1.4 Mq. 

In our work, we use the initial masses of NSs via CCSN in ? and M^ s = IAMq in different 
simulations, respectively. For NSs via AIC, following ?, we take M^ s = 1.3M0. Similarly, for NSs 
via EIC, we also take M^ s = 1.3MQ. 

In addition, nascent NS receives additional velocity ("kick") due to some still unclear process 
that disrupts spherical symmetry during the collapse or later Dichotomous nature of kicks which 
was suggested quite early by ?. Observationally, the kick is not well constrained due to numerous 
selection effects. Currently, high kicks (~ 100 km s" 1 ) are associated with NS originating from 
CCSN, while low kicks (~ 10km s -1 ) with NS born in EIC and AIC (?). We apply to core-collapse 
NS Maxwellian distribution of kick velocity 

PK) = #4^ /2 < (i) 

Variation of o"k between 50 and 200 km s -1 , introduces an uncertainty < 3 in the birthrate of low- 
and intermediate-mass X-ray binaries (?). ? discussed the effects of parameter on LMXBs' 
populations. Since in this paper we focus on the physical parameters that mostly affect the masses 
of NS, we do not discuss the effects of on SSs' population. Following ?, we take = 190 km 
s- 1 in CCSN, and a k = 20 km s" 1 in EIC and AIC. 

2.2. Mass of Accreting Neutron Star 

In LMXBs, NSs accrete the matter from their companions via Roche lobe flows or stellar winds. 
The interaction of a rotating magnetized NS (single or in a binary system) with surrounding matter 
has been studies by many authors (?????, e. g.,). 
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Using a convenient way of describing NS evolution elaborated by ? and a recent model for 
quasi-spherical accretion including subsonic settling proposed by ?, ? gave a detailed simulations 
for spin period evolution and matter accretion of NSs in binaries. In this work, we adopt their 
model. ? assumed that all matter transferred is accreted by the NS in an accretor stage. We 
introduce a parameter /3 which is the fraction of transferred matter accreted by the NS, and the 
rest of the transferred mass is lost from binary system. The lost matter takes away the specific 
angular momentum of the prospective donor. The value of j3 has been usually set to 0.5 (???). In 
our work, we set (3 = 0.1, 0.5 and 1.0 in different simulations. 

2.3. LMXBs with Nascent SSs 

As noted in when the gravitational mass of a NS reaches to ~ 1.8Mq, the NS can turn into 
a SS via the hydrodynamically unstable conversion or the slow conversion. In the former mode, 
the binary system may been disrupted (?), and it becomes two isolated star. One of them is an 
isolated SS. However, it is difficult to know the effects of the hydrodynamically unstable conversion 
on binary systems. Therefore, in our work, we consider two extreme cases: (i) In order to simulate 
all potential properties of qLMXBs, we assumed all LMXBs are not affected and survive after the 
conversion, that is, LMXBs become qLMXBs when the masses of accreting NSs are larger than 
1.8Mq| (ii) We assumed that all LMXBs are disrupted after the conversion, that is, there is no 
qLMXBs. However, we can discuss the origin of isolated submillisecond pulsars. 

In our simulations, if a nascent NS has larger mass than 1.8 M@ it is a SS. Therefore, in the 
paper, qLMXBs include some LMXBs in which the nascent NSs have larger mass than 1.8 M . ? 
suggested that SSs could have high maximum masses (See Figure 2 in ?) if they are composed of 
the Lennard- Jones matter. In our work, we assume that the maximum mass of SS is 3.0 M@. 

3. Results 

We use Monte Carlo method to simulate the initial binaries. For initial mass function, mass- 
ratios, and separations of components in binary systems, we adopt the distributions used by us in 
??. We assume that all binaries have initially circular orbits. After a supernova, new parameters of 
the orbit are derived using standard formulae, (e. g., ?). It is well known that theoretical models 
of the population of LMXBs depend on badly known input parameters, such as kick velocity and 
common envelope treatment (e. g., ??). However, in this pioneering study of qLMXBs we focus on 
the effects which are important for the masses of NSs: the efficiency accreted of transferred matter, 
(3 ( (3 = 0.1, 0.5 and 1.0 ), and the mass of nascent NS via CCSN (See Figured] ). We use 1 x 10 8 
binary systems in our Monte-Carlo simulations. 
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3.1. Mass Increases of Accreting NSs 

Figure [2] gives the mass increases of accreting NSs in LMXBs. According the assumptions of 
SSs formatting from NSs with masses higher than 1.8 M , l%o (/3 = 0.1) — 10% (/3 = 1.0) of 
LMXBs are qLMXBs in our simulations. This proportion greatly depends on input parameters /3 
and the initial masses of nascent NSs. In the cases of (3 = 1.0 and /3 = 0.5, most of SSs in qLMXBs 
come from low-mass NSs with 1.4 M@, that is, they have accreted ~ O.4M0 matter. In the case of 
f3 = 0.1, most of SSs in qLMXBs originate from CCSN. 

Mass increases of accreting NSs in LMXBs depend on not only input parameter /3 in this work, 
but also orbital periods and NSs' companions. ? showed that most of NSs in LMXBs with WD 
donors have low mass- accretion rates (~ 10 _12 MQyr~ 1 ) and most of LMXBs with WD donors are 
transient. Therefore, the masses of accreting NSs in LMXBs with WD donors hardly reach to 1.8 
M . Less than 10% (in case of /3 = 0.1)— l%o (in case of /3 = 0.5, M^ s = 1.4M ) of qLMXBs have 
WD donors in our simulations. Most of qLMXBs have MS donors. From now on, we just discuss 
LMXBs and qLMXBs with MS donors. 

3.2. Properties of qLMXBs 

According to our assumption that the conversion of NS matter to SS matter does not affect 
binary systems, we can simulate some observational properties of LMXBs and qLMXBs, and wish 
to find some differences between them. Then, taking the case of /3 = 0.5 for example, we discuss 
the some properties of qLMXBs. 

X-ray luminosities (mass-accretion rates), spin periods and orbital periods are important pa- 
rameters of LMXBs. Figure [3] shows the accretion rates by NSs or SSs in LMXBs and qLMXBs, 
and the X-ray luminosities which are approximated as 

L x = t/Mnsc 2 ~ 5.7 x 10 35 erg s" 1 ^) x U-j^-i ), (2) 

where rj ~ 0.1 is the efficiency of converting accreted mass into X-ray photons. We can find that 
there is not significant difference between accretion rates by NSs or SSs in LMXBs and qLMXBs. 

The spin periods of the rotating magnetized NSs mainly depend on their mass-accretion rates 
(?). If accreting SSs are similar to NSs, we can simulate the spin periods of accreting SSs. Figure 
U] gives the distribution of the spin periods of NSs in LMXBs or SSs in qLMXBs. In general, spin 
periods of SSs in qLMXBs are shorter than those of NSs in LMXBs. 

Figure [5] shows the the distributions of the orbital periods P Q rb of qLMXBs and LMXBs vs. 
masses of their secondary stars. Similarly, there is not significant difference between qLMXBs and 
LMXBs. 

As the above descriptions show, it is very difficult in our model to differ from LMXBs and 
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Fig. 1. — Masses of nascent NSs via CCSN vs. stellar initial masses. The solid line comes from ?, 
and the dashed line means that masses of nascent NSs are equal to 1.4M . 
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Fig. 3. — Distributions of accretion rates (X-ray luminosities) by NSs or SSs in LMXBs and 
qLMXBs, respectively. 









1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 








qLMXBs 








LMXBs 








_ 

1 

— i 
1 








L _ 

L _ 



Fig. 4. — Number distributions of the spin periods of NSs or SSs in LMXBs and qLMXBs, respec- 
tively. The numbers are normalized to 1. 
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qLMXBs. The most effective way is to measure the masses of compact objects in LMXBs if the 
assumption that NSs with masses larger than 1.8 Mq are SSs is right. 

? suggested that compact object in 2S 0921-630 (It is a LMXB) has a mass between ~1.9 — 2.9 
Mq. According to our assumption, the compact object is a SS. However, ? considered that ? 
overestimated the rotational broadening and the mass of compact object in 2S 0921-630 is ~1.44 
M Q . The orbital period of 2S 0921-630 is 9.006±0.007 days, and its X-ray luminosity is ~ 10 36 erg 
s _1 (?). Results of simulating LMXBs and qLMXBs in our work both cover these observations. 
Therefore, we can not conclude whether compact object in 2S 0921-630 is NSs or SSs. 

? gave that PSR 1614-2230 has a mass of 1.97 Mq. And, it is a millisecond radio pulsar 
(Pulsar spin period is 3.15 ms) in an 8.7 day orbit, and its companion has a mass of 0.5 M a dot. 
Although PSR 1614-2230 is not X-ray binary, it may come from X-ray binary. ? suggested that 
PSR 1614-2230 descended from an LMXB very much like Cyg X-2 (P orb =9.8 days, M NS = 1.7M 
and M2 = 0.6M Q , see ?). Figure [5] shows the distributions of the orbital periods -P or b of qLMXBs 
and LMXBs vs. masses of their secondary stars. Our results cover the positions of orbital periods 
and companion masses of Cyg X-2 and radio millisecond pulsar binary PSR 1614-2230. Our work 
support that PSR 1614-2230 originate from a LMXB. Now, PSR 1614-2230 may be a SS. 



3.3. Submillisecond Pulsars 

? suggested that an isolated submillisecond pulsar spinning at ~ 0.5 ms could strongly hint 
the existence of SSs. If the conversion of NS matter to SS matter is hydro dynamically unstable and 
LMXBs are disrupted, the nascent SSs are isolated. Assuming that one binary with M\ > O.8M0 is 
formed per year in the Galaxy (???), we can estimate that the occurrence rate of hydrodynamically 
unstable conversion is about 5 — 70 per Myr. 

As Figure [6] shows, the majority of NSs at the beginning of the conversion have spin periods 
longer than 1 ms. If the angular momentum is conserved and no unknown physical mechanic 
spins up NSs and nascent SSs during the conversion, the spin periods of nascent SSs depend on 
the change of moment of inertia. ? considered that a typical ejected mass of hydrodynamically 
unstable conversion is ~ 1O _3 M0. Therefore, the change of moment of inertia is determined by 
the difference of radius between NS and SS. ? estimated that NS could shrink by as much as 30%. 
Then, many nascent SSs have spin periods of ~ 0.5 ms in our simulations. However, this result 
greatly depends on the equation of state of NSs and SSs which are poorly known. If NS only shrinks 
by as much as 10% (Private discussion with Xu), it is difficult for the nascent SSs to spin up to 0.5 
ms. 
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Fig. 5. — Similar to Figured but for the distributions of the orbital periods P or b of qLMXBs and 
LMXBs vs. masses of their secondary stars, respectively. 





Fig. 6. — Number distribution of the spin periods of NSs at the beginning of the conversion of NS 
matter to SS matter. 
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4. Conclusions 

Employing the population synthesis approach to the evolution of binaries and using the inter- 
acting model of a rotating magnetized NS with surrounding matter, we investigate the mass change 
of NSs in LMXBs and the possibility from NSs converting SSs in LMXBs. Our results shows that 
about l%o — 10% of LMXBs can produce SSs. These SSs may exist in qLMXBs or be isolated, 
which depends on physical model of the conversion of NS matter to SS mater. 

Our toy model can not conclude whether there are SSs in the Galaxy and can not give what 
properties qLMXBs have. In further work, we need detailed physical model (equation of state 
about NS and SS and the conversion of NS matter to SS mater) to improve our work. 
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